Abstract 



Considering the spin-coefficient version of the ieft-fiat vacuum Einstein 
equations, aii but one of the fifty equations can be explicitly integrated 
via the introduction of five spin-weight s=-2 complex potentials. The 
final equation is a non-linear wave equation for the last of the potentials. 
Solutions to this equation determine solutions for the entire system. 

Solutions for several special cases are obtained. 



1 I. Introduction 

1.1 Background 

The self-dual/left flat (or anti-self dual/right /flat) Ricci-flat equations and their 
spaces with associated metrics has been a subject of considerable interest since 
the late 1970s. The very large literature has been associated with a variety of 
linked interests. 

They play a fundamental role in the properties of real asymptotically flat 
space-times via solutions to the 'good-cut' equation 1,2,3,4,5,6,7 , they define Pen- 
rose's asymptotic twistor space 8,9,10 , they (the solutions) are intimately related 
to (in fact are) Penrose's non-linear graviton fields, the equations themselves 
are studied as a beautiful example of a non-linear integrable system 11,12,13 and 
finally the Euclidean versions, the so-called gravitational instantons 14,15 (the 
general relativistic generalization of the Yang-Mills instantons) have played a 
role in attempts at understanding quantum gravity. 

From the beginning there have been a variety of attempts at integrating the 
field equations and even now - 40 years after their popularity began - only a 
relatively small number of solutions appear to be known. The first solutions - 
a very limited number - were obtained via the 'good cut' equation. The best 
know is the Sparling- Tod metric 6 / Eguchi-Hansen 7 metric. Many others were 
discovered by imposing special conditions as for example algebraic specialness 16 
or (Killing) symmetries 17 . In addition there is a large literature on self-dual 
gravitational instantons 14 . 

To our knowledge there have been very few attacks on the problem of finding 
(or even of studying) the general solution of the self-dual Einstein equations on 
general complex four dimensional manifolds - either globally or locally 4 . It is 
the purpose of this note to return to this issue. 

Although all H-spaces (i.e., spaces arising from solutions of the good-cut 
equation) are self-dual spaces, it is not yet clear whether all self-dual spaces are 
H-spaces. Nevertheless we will not make a distinction here between the two. 



1.2 Modus Operandi 

We begin, in Sec. II, with the fully general Einstein equations in the spin- 
coefficient formalism written as complex equations on a thickened region of 
C 4 in a neighborhood of the real 'slice' R . In the standard "real" version of 
the formalism, the basic field variables, i.e., the spin-coefficients, the metric 
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variables and Weyl tensor components, are, in general, complex. In the field 
equations themselves, their complex conjugates appear symmetrically so that 
the equations themselves, as a complete set are real. In our present (SD) ver- 
sion however, the previous "complex conjugates" are now completely freed up 
from their conjugate counterparts and are independent variables. These field 
equations are then reduced by setting the anti-self-dual part of the Weyl tensor 
to zero, thereby imposing the self-dual condition on the equations. 

A coordinate system is introduced, in Sec. II, by choosing a world-line, £, 
in the thickened region of C 4 : with the (complex) u, being the 'complex time' 
along the world line. The null cones with apex on the line are labeled by 
u. The complex sphere of null directions at the apex of each cone (and their 
associated null geodesies) are labeled by the complex stereographic coordinates 
(C, C) while the affinc parameter along the null geodesies is given by r. Note 
that both u and r have values close to the real while £ is close to the complex 
conjugate (. 

From these imposed conditions we easily show, Sec. II, that the complex 
divergence, p, of the null geodesies with apex on £ is given by p = — r _1 and 
the left shear by a = 0. 

At this point the entire set of field equations (about 50) are divided into 
three sets: (1) the radial equations which involve the r derivative, (2) those 
involving the angular derivatives, i.e., (C, C) an d (3) those that contain the 
time-derivative, u - about (14) of them. 

1.3 Results 

What is a bit surprising is that every one of the r-derivative equations 
can be exactly integrated in terms of four spin- weight (s = —2), potentials, 
(7oi 7ii 72j 73 7 74)7 where each can be expressed as r-derivatives of the next one, 
so that the only independent one is j 4 . 

The angular equations then are then used to establish relationships 
between some of the free "constants" of the r- integrations. 

Finally when all these results are inserted into the (14) time-derivative 
equations we discover that they are all closely related and can be reduced to 
a single complex spin-weight ( s = —2), non-linear wave equation for the last 
potential 74. This last equation carries all the evolution information for the 
entire set of equations. 

Though we do not know of a way to give general solutions to this equation, 
it is quite easy to produce many special solutions. 

2 II. The Field Equations 

2.1 Spin-Coefficient version of Einstein equations 

The familiar spin-coefficient version of the Einstein equations 18 - with their "con- 
jugates" explicitly inserted and written out in detail - are: 
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The metric equations: 



Al a -Dn a = (j+^)l a + (e + e)n a -(T+W)m a -(T + TT)m a , (1) 

Sl a -Dm a = {a + f3-Tf)r + Kn a -crm a -(]5 + 6-e)m a , 

5l a ~Dm a = (a + P-7r)l a +Kn a -Wm a -(p + e-e)m a , 

5n a -Am a = -vl a + (r -75 - /3)n a + Xm a + (p - 7 + j)m a , 

5n a -Am a = -Vl a + (T-a-~j3)n a + \m a + (-p-^ + -/)M a 

~8m a - 8m a = {p - p)l a + (p - p)n a - (75 - (3)m a + (a - ~p)m a . 

The spin-coefficient equations: 

AX-Su = -(p + 7J+37-7)A+(3a + ^ + 7r-T)z/-*4 (2) 

AX-5V = -(p + 71+ 37-7)!+ (3a + (3 +W - t)V - * 4 

dp — 5a = p(75 + (3) — a(3a — /?) + (p — p)r + (p — p)fc — "fi 

5p - 6a = p(a + /?) - ct(375 - /?) + (p - p)r + (p - p)7c - *i 

5a -J/3 = pp - Act + aa + /3/3 - 2a(3 + ^{p - p) + e(p - p) - # 2 

(575 — 5(3 = JTp — Act + aa + — 275/3 — 7(p — p) — e(p — p) — ^2 

<5A-Ip = (p - p> + (p - 7J)7r + p(a + /3) + A(75- 3/3) - * 3 

~SX-5ji = -(p- p)z7- (p - 7J)7f + p(75 + /3) + A(a - 3/3) - * 3 

(5i^ — Ap = p 2 + AA + p(7 + 7) — Vn + v(r — 3(3 — 75) 

5u — Ap = ~p 2 + XX + 71(7 + 7) — vlx + 77(r — 3/3 — a) 

57 — A/3 = 7(r — 75 — /3) + pr — av — eV — /3(7 — 7 — p) + aX 

57 — A/3 = 7(t — a — (3)+ JFF — av — ev — (3{j — 7 — p) + 75A 

5t — Act = pa + pX + r(r + (3 — a) — a(3j — 7) — nv 

St — Aa = pa + pA + r(r + /3 — a) — 7f (37 — 7) — ~nv 

Ap — 5t = —pJJ. — ctA + t{(3 — a — t) + (7 + 7)p + kv — ^2 

Ap — <5t = — pp — 7f A + r(/3 — 75 — t) + (7 + 7")p + 7577 — ^2 

Aa-Sj = v(p + e) - A(r + /?) + 0(7 - 77) + 7^ - r) - * 3 

Aa - 5j = V(p + e) - A(r + /3) + 75( 7 - p) + t"(/3 - r) - * 3 
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Dp- 


- OK 




p 2 + (7cr + (e + e)p — kt — K(3a + j3 — tt) 




Dp- 


r — 

- OK 




p 2 + (7CT + 


(e + e)p — kt — K(3a + j3 — n) 




Da - 


- OK 




(p + p)(7 + 


(3e - e)(7 - (r - W + a + 3/3)k + * 




Da - 


T — 

- OK 






(3e - e)a - (r - tt + a 4- 3/3)7c 4- *o 


(3) 


Dr — 






(r + lf)p + (r + 7r)(7 + (e - e)r - (37 + 7)k + *i 




Dr — 


A — 




(r + tt)p + (t + tt)ct - (e - e)r - (37 + j)k + *i 


(4) 


Da - 


T 

- oe 




(p + e- 2e)a + pa - ~/3e - kX - «7 + (e + p)tt 




Da 


c— 

- 06 




(p + e — 2e)a + /3<7 — pe — kX — kj + (e + p)7r 


(5) 


Df3 - 


- 5f 

L/C 




(a + ir)a 4 


- (p - e)/3 - (/x 4- 7)k - (a - 7f)e 4- *i 




DP 


-5e 




(a + W )a 4 


- (p - e)/3 - (m 4- 7)k - (a - 7r)e + *i 


(6) 


£>7 - 


Ae 




(r + 7f)a 4 


- (r 4- 7r)/3 - (e + e)7 - (7 + 7)e 4- T7r - 


W+ *2 


L>7 - 


Ae 




(r + 7r)a + (t + 7r)/3 -(e + e)7 - (7 + 7)e 4- T7r - 


+ * 2 (7) 


DX - 


- Sir 




pX + ap + 


tt 2 + (a - /3)tt - vK - (3e - e)A 




DX- 


- 5lf 




pX + <T~p + 


7f 2 + (a — /?)7f — Fre — (3e — e)A 


(8) 


Dp- 


- Sir 




pp + aX + 


7T7T — (f +?)/!- 7T ((5 — /3) — I^K 4" *2 




Dp- 


- Sir 




pp + aX + 


7T7f — (e 4- e)p — 7f (a — f3) — vk + *2 


(9) 


Dv- 


Att 




(r 4- 7r)p + (t 4- 7f)A + (7 - 7)tt - (3e + e)v + * 3 




DV- 


Att 




(r + 7f)7i + (? + tt)A - (7 - 7)tt - (3? + e)z7 + * 3 


(1UJ 


and finally the Bianchi Identities: 








(5*o 


- £)*i = 


(4a - 7t)* - 2(2p 4- e)*i 4- 3k* 2 , 


(11) 






5*! 


-£>* 2 = 


A* 4- 2(a - 7r)*i - 3p* 2 4- 2k* 3 , 








5*2 


-£>*s = 


2A*i - 3tt* 2 4- 2(e - p)* 3 4- k*4, 








5*3 


- Z)*4 = 


3A* 2 - 2 (a 4- 2tt)* 3 4- (4e - p)* 4 , 








5* 


- £>*i = 


(4a - 7f)* - 2(2p 4- e)*i 4- 3k* 2 , 


(12) 






5*1 


- £>*2 = 


A* + 2 (a - 7f)*i - 3p*2 4- 2k* 3 , 








5*2 


-£*3 = 


2A*i - 37f* 2 + 2(e - p)* 3 4- k* 4 , 








5*3 


- Z?*4 = 


3A* 2 -2(5+ 27f)* 3 4- (4e - p)* 4 , 








A* 


-5*! = 


(47 - p)*o - 2(2t + /3)*i 4- 3cr* 2 , 


(13) 






A*i 


- 5*2 = 


+ 2(7 - p)*i - 3r* 2 4- 2cr*3, 








A* 2 


- 5*3 = 


2^*! - 3p* 2 4- 2(/3 - t)* 3 4- (7*4, 








A* 3 


-5*4 = 


3^*2 - 2(7 4- 2p)* 3 4- (4/3 - r)* 4 . 
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AWo-JWi = (47-7Z)* -2(2r + /3)*i +3^*2, (14) 

A*i-M 2 = + 2(7 -7Z)*i - 3r* 2 + 20^3, 

A* 2 -J*3 = 217*! - 37**2 + 2(/3 - r)f 3 + 0^4, 

A* 3 -M 4 = 3F* 2 - 2(7 + 2/Z)* 3 + (4/3- r)* 4 . 

The A" = (l a ,n a 7 m a 7 m a ) are the components of a null tetrad, the * 's 
and * 's are the components of the Weyl tensor, A"V a = (D, A, 5, 5) are the 
directional derivatives. All the other variables are the spin-coefficients. The 
Einstein equations are already built into the system by virtue of the Ricci tensor 
having been set to zero. 



3 III. Left-Flat Spin- Coefficient Version of Ein- 
stein equations 

There are now several different conditions that are now imposed on these equa- 
tions to restrict them to left fiat equations and in addition to greatly simplify 
them. 

The first - and most basic - is to impose the left-fiat (or self-dual) condition 
on the equations. This is accomplished by simply taking 

(* ,*1,*2,*3,*4) = 0, (15) 

which is equivalent to imposing on the Weyl tensor C a bcd that 

C abcd e cd *t = iC ab ef . (16) 

All the barred variables (e.g., a, A....) are independent and freed from their 
dual counterparts. This is denoted by replacing the bars by tildes, (e.g., a, A...). 

The other conditions (both coordinate and tetrad conditions) are imposed 
on the tetrad vectors, A". We choose a world- line, I, in the complex thickened 
i? 4 and a one parameter family of (complex) light cones with apex on I labeled by 
the complex coordinate x° = u. The null generators of the cones (null geodesies) 
are labeled by the points of the thickened sphere with complex stereographic 
coordinates x A — (£, £), with C ~ C The affine parameter of the null geodesies 
is the complex radial coordinate x 1 — r. These conditions allow us 18 to write 
the directional derivatives (or tetrad) as 

(17) 
(18) 
(19) 
(20) 



D = 


dx a 


d 
dr 


V = 


"•if- 

dx a 




5 = 


dx a 


* 9 + c 4 d 

dr dx A ' 


5 = 


dx a 


„ d + ~A d 
dr dx A ' 
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and associated metric, [g ab = l a n b + n a l b — m a m b - m a m b ], 

(21) 
with 

2{U-uw), (22) 
X A -(wt A +u,t), 

By appropriate tetrad transformations 18 , (parallel propagation of n a , m a , m a 
and appropriate scaling of l a ) we can put the following conditions on the spin- 
cocmcients: 

k = it = e = k = ir =1= = p — p = 0. (23) 
We also have, but do not use immediately, 

r = 5 + /3, 7 = a + /3. (24) 

These conditions are then inserted into our Einstein field equations to obtain 
our final set of field equations. The equations are displayed in groups: first, (17) 
of them, that contain the D derivatives, then, (7), containing only the angular 
derivatives, {5,5 ) and finally, (20), containing the time derivative, A. 



g ab = 



o , 1 , 

1, 9 U : 

0, g 1A 





9 1A 



.,22 



9 2A = 



6 



Dp 




p 2 + aa 


Da 




2pa, 


Da 




2pa + *o 


Dt 




rp + fa, 


Dt 




Tp + Ta + *i 


Da 




pa + 


Da 




p5 + /3<7 


DB 




ctcr -+- pB 


Dp 




aa + p/3 + *\ 


D7 




ra + r/3 


Dj 




t5 + t/3 + $ 2 


DX 




pA + aji 


DX 




pX + ajl 


Dii 




pp, + aX 


Dp, 




pp + aX + * 2 


Dv 




rp + tA 


DV 




rp + ?A + $ 3 



DU 

DX A 
Duj 
DZj 

D^ A 

~A 

Dt; 



TUJ + TLJ — (7 + 7) , 
~A 



= au> + poj — t, 

= auj + puj — r, 
~A A 

~ A ~ A 



(5* - D^>i = 45*o -4p*i, 

(5*i - L>* 2 = A* + 25*i - 3p* 

(5* 2 - £>* 3 = 2A*i - 2p* 3 , 

(5*3 - L>* 4 = 3A* 2 - 25*3 - P* 



Sp-Sa = p(a + /?) - cr(3a - /?) (49) 

5p-<55 = p(a + /3) -5(35-/3) -*i (50) 

(5a — S(3 = pp — Act + a5 + /3/3 — 2a(3 (51) 

55 - sji = pp - A5 + aa + /3/3 - 25/3 - * 2 (52) 

<5A-dp = p{a + /3) + A(5 - 3/3) (53) 

5A — 5^1 = p(5 + /3) + A(a - 3/3) - * 3 (54) 



8uj — Slu = (fj, — p) — (a — P)uj + (a — f3)u>. (55) 

S£ A -sf = -{a-P)t + {a-~^ A (56) 

SU - Auj = -u + Xuj + (p- 7 + 7)0;, (57) 

SU - Aw = —v + Xuj + (p — 7 + 7)w, (58) 

M^-Af 4 = Af + ( M -7 + 7)^ (59) 

5X A -A^ = K A + {p-l + l)T (60) 



AA 


-Sv = 


-(/i + p)A-(3 7 -7)+2aA 


(61) 


AA 


-Sv = 


-(p + p)A - (37 - 7) + 25A - * 4 


(62) 


Sv- 


- Ap = 


p 2 + AA + p( 7 + 7) - 2/3^ 


(63) 


~8v- 


- Ap = 


p 2 + AA + p( 7 + 7) - 2/3£ 


(64) 


5 j - 


-A/3 = 


pT — av — /3(7 — 7 — p) + aA 


(65) 


<5 17 - 


-A/3 = 


pr — av — /3(7 — 7 — /I) + 5A 


(66) 


<5r- 


-Act = 


(/i(j + pX) + t(t + (3 — 5) — er(37 - 7) 


(67) 


55- 


-A5 = 


(p5 + pA) + r(r + /3 - a) - 5(37 - 7) 


(68) 


Ap 


-St = 


-(pp + o-A) + t(/3 - a - t) + (7 + 7)p 


(69) 


Ap 


-St = 


-(pp + 5A) + 5(/3 - 5 - r) + (7 + t> - * 2 


(70) 


Aa 


— 8 j = 


^p - A(r + /3) + a (7 - p) + 7(/3 - 5) 


(71) 


A5 


— 8j = 


vp - A(r + /3) + 5(7 - p) + 7(/3 - r) - $3 


(72) 



A* -5*i = (47 - p)$ - 2(27 + /3)$!+ 35*2, (73) 

A*i-5* 2 = + 2(7 -p)$i - 3r* 2 + 25*3, (74) 

A* 2 -5* 3 = 2?*! - 3p$ 2 + 2(/3 -5)$ 3 + 5*4, (75) 

A* 3 -5* 4 = 3£* 2 - 2(7 + 2p)* 3 + (4/3 -r)* 4 . (76) 
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In the following section we integrate explicitly all the D equations. The 
solution to each equation will have an associated 'constant' (actually a function 
independent of r) of integration. These 'constants' are determined by either the 
angular equations or by the fact of all the null surfaces being light-cones with 
their origin on a world- line. The evolution equations (i.e., those with A) turn 
out - initially - to be extremely complicated to deal with. Nevertheless - and 
quite surprising - in the end they all are equivalent to one single non-linear wave 
equation. 



4 IV. Integration 

4.1 Radial Integration 
4.1.1 preliminary integration 

The first thing to note is that Eqs. ([25"T) and (|26|) can be integrated immediately, 
(with the coordinate origin for r taken at the apex of the light-cone) as 

P = --■ (77) 
r 

a = 0. (78) 

The later follows from the fact that a vanishes at a light-cone apex. 

The field equations are greatly simplified by these results. The spin-coefficient 
Eas. ([28l3Tl l32l 

|3"T|) I3"51 decouple from the remainder and using the results from fiat space light- 
cones with origins on a world-line 1 we have 



T = 


a + /3 = 0, 


(79) 


a — 




(80) 


P = 


—a, 


(81) 


A = 


o, 


(82) 


P = 


-r- 1 . 


(83) 



The form of the function P, (I80|) , which determines the 2-sphere metric at 
the cones apex via Eqs.(|22p. is set as a coordinate condition. 

In a similar fashion, with r = 0, the metric Eqs.([44j>. ((46]) integrate to 

t 4 = r- 1 (e C ,e~ C ) = -r- 1 (P,0) > (84) 
lj = 0. (85) 
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4.1.2 The Radial Bianchi Identities & Weyl Tensor 



Turning now to our first non-trivial integrations, using the preceding results, 
the four radial Bianchi Identities, (|48a[) - l48df become 



d r *l = 


-4r 






d r *2 = 


-3r 




-r _1 3#i 


a r $ 3 = 


-2r 






5 r $ 4 = 


—r~ 


^4 ~ 





(86a) 
(86b) 
(86c) 
(86d) 

where Eqs. ([T7j) and (|20l) . along with 3?7 S = P 1 ^ s -§^(P s r] s ), were used. 
Writing these equations, with n — 1,2, 3, 4, succinctly as 

9r*5-n = -nr- x * 5 -n - r- x 8*4-n (87) 
the solutions are given by 

§ 5 -„ = ^ n $5°-n - r _n 3 / r"- 1 $ 4 _„dr, (88) 

the integral being indefinite, with ^^l n the " constants of integration" . 

Our integration process continues with the introduction of the set of five, 
spin-weight s = —2, potentials, T5_„(r), (n = 1,2,3,4,5) by 

*6-n = (-l)"- 1 r-"+ 2 g( 5 -")9 r T 6 _„ (89) 



and inserting them into the integral of Eq. 
Integration leads to 

$ 5 _„ = r-"*5 2 n(0) - (-l)"r-"9( 5 -")T 5 _ n (r), (90) 
where, from regularity at r = 0, we have 

*S%) = °- ( Q1 ) 

Remark 1 W^e emphasize that there are two different types of situations. In 
one case the potentials, T, vanish sufficiently fast at r = so that the origin is a 
regular point of the manifold. In the other case there are intrinsic singularities 
hidden in the T(r) at r = 0. See the following section for examples.] 

Written explicitly, Eqs. (|9"0)) are 

(92) 
(93) 
(94) 
(95) 
(96) 



*0 


= r 


" 5 T , 


*1 




•r _4 9Ti 




= r 


~ 3 g 2 t 2 , 


$ 3 






*4 


= r 


- 1 g 4 f 4 . 
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By equating Eqs. ([8"9"|) and (|M|) . we find the 'ladder' relationship between the 
different potentials, 

T 5 _ n - r 2 a r T 6 _„ (97) 

or 

T = r 2 a r T 1; (98) 

Ti = r 2 d r t 2 , (99) 

T 2 = r 2 a r T 3 , (100) 

T 3 = r 2 d r t 4 . (101) 

Knowledge of T4 (which we can consider as free initial data) determines 
the others or equivalently, from To with the four constants of integration, from 
going up the ladder to T4, determines the others. 

In the case of the regularity of the Weyl components at r — 0, we must 
impose conditions on the T„ so that near r = 

T n = 0(r 5 ~"). (102) 

This, with suitable differentiability, follows from T4 = Ar + 0(r 2 ). 

For conventional asymptotic flatness (peeling) we require that, as r — ► 00, all 
the potentials tend to a constant, T„ = c n + 0(r~ r ). 

From Eq. (|9"7} we have the asymptotic behavior of the potentials: 



To 


= r 2 a r T 1=Co 


+ 0(r- 




Ti 


= r 2 a r f 2 = ci 


- c r~ 




T 2 


= r 2 d r T 3 = c 2 


— C\T~ 




T 3 


= r 2 <9 r T 4 = c 3 


- c 2 r~ 




T 4 


= c 4 - c^ 1 + 


1 

2 C2 ^ 


-2 _ 



(103) 

0(r- 2 ), (104) 
ic o r- 2 +0(r- 3 ), (105) 

\cir~ 2 -jUr" 3 + 0(r- 4 ), (106) 

i Cl r- 3 + ^ C or- 4 + 0(r- 5 ). (107) 



4.1.3 Remaining Radial Equations 

The remaining spin-coefficient and metric equations can be integrated and ex- 
pressed in terms of the potentials. 

We first investigate Eq.(|27|) using v&o = r 5 ^o, i-e., 



Do = -2r- 1 ? + r- 5 T . (108) 

Its solution is given by 

a = a°r- 2 + r~ 2 [ r _3 Todr. (109) 



with (again) an indefinite integral and ct° being the constant of integration. 
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The integral term can be greatly simplified via repeated use of Eq. ipT)) . 
From 

J r- 3 t dr 



we have 

a = a°r- 2 +r- 3 t 1 +r- 2 t 2 (111) 

From the vanishing of the T at r = and the regularity of the light-cones 
at the origin, we take a = 0, so that 

5 = r~ 3 T 1 +r- 2 T 2 . (112) 

Using the same methods, i.e., the formal radial integration of the equations 
followed by the repeated use of the ladder relations, Eq. (l97l) . all the remaining 
radial equations can be integrated and expressed in terms of the T„. Many of 
the associated "constants of integration" are determined by the behavior near 
r = 0, several are determined by use of the angular equation. (It is likely that 
even those determined from the angular equations could have been determined 
from their r — behavior.) 

The following is the full set of solutions to all the spin-coefficient and metric 
radial equations. 



r -1 r fidr = J [0 r (r -1 fi) + r- 2 t 1 ]dr (110) 

/r pr 
r~ 2 t 1 dr^r- 1 T 1 + I d r t 2 dr 

r-ifj +T 2 
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Spin-coefficients: 



a = T = a + /3 = 0, r = a+ $, (113) 

P = P = -r-\ (114) 

o = r- 2 T 2 + r- 3 T 1; (115) 

t = -r- 2 dt 2 -r- 1 dt 3 , (116) 

a = r^a = --r~ 1 d i P, (117) 

(3 = r- 1 f3° = -a, (118) 

a = -r-i^P-dtPi^+r- 2 ±d c Pf 2 , (119) 

= r-^^P-Bts-dcPt^-r- 2 (3T 2 + ia c PT 2 ) (120) 

7 = -a c pgT 4 -r- 1 ia c pgT 3 (121) 

7 = g 2 T 4 + 9c^T4 + r- 1 (g 2 T 3 + ia c pgT 3 ) (122) 

A = (123) 

fj, = -r- 1 (124) 

A = - 2r~ 1 T 3 -r- 2 T 2 (125) 

Jt = -r- 1 +r- 1 3 2 f 3 (126) 

v = -9 3 T 4 (127) 

v = 2dt i + r- 1 dt 3 (128) 

Metric Variables: 

Z A = (£ C ,&=r- i e A = -r- 1 (P,0), P=1 + CC (129) 

t = (f ,? C ) = -r- 1 (0,P)-(P,0)(2r- 1 T 3 + r- 2 T 2 ) (130) 

w = r _1 w° = (131) 

w = r" 1 ^ + 3f 3 = 9T 3 (132) 

AT A = (X c ,^) = (P,0)[r _1 3t3 + 29f4] (133) 

C/ = -l-rg 2 T 4 (134) 

From these results we have the derivative operators which are needed in the 
angular and evolution equations, namely 

D = d r , (135) 

V = 9„-(l + rg 2 f 4 )5 r + (r- 1 Pgf 3 + 2Pgf 4 )5 c , (136) 

S = -r _1 P0 c , (137) 

5 = df 3 d r -r- 1 Pd^-P(r- 2 t 2 + 2r- 1 T 2 )d < . (138) 
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4.2 Angular Equations 

Virtually all the angular derivative spin-coefficient equations turn out to be 
identities when the results of the previous section are used. Some determine 
the "constants of integration" but they are almost certainly also determined 
by the behavior near r — 0. Several examples, (already used in the previous 
subsection) are: 

£V + 7-7) + Ae A =>7° = 0, (139) 
Cp~a)t A + {a-P)f =>a° = -^d c P, (140) 

- a)uj + (a - (3)uj + (fj, - fx) => Identity, (141) 
+ 1 + \u -D => p° = 0. (142) 

4.3 Evolution Equations 

Analyzing the twenty evolution equations for the five potentials, Eqs. (|5Tj) - (fTo| , 
initially presented a very difficult challenge. The process of inserting all the 
previous results, i.e., Eqs. (j92j) - (j96]) and Eqs. (I113|) - (I134|) . into the evolution equa- 
tions was a daunting task. Analyzing any one could take several days. The 
algebra was long and complicated with the easy production of many errors. 
The saving aspect of it was the fact that the equations were interrelated, in the 
sense that final forms of many of the equations could be compared or reduced 
to others, thereby giving a method of checking consistency and thus finding the 
errors. 

Towards the end of the analysis all the evolution equations depended on just 
the four evolutionary Bianchi Identities, Eqs.(f57]) - (jM)) . They however, using 
Eqs. (|98|) - (|101|) . are further reduced to a single equation for the last of the 
potentials, i.e., T 4 . 

Taking the four evolutionary Bianchi Identities, Eqs. (j73|) - (|76|) . substituting 
the results of the radial integrations - with lengthy calculations using Eqs. (f9"5]) 
- (jlOip . - we obtain the following results: 

To - 93T! - 2T 3 g 2 Ti + 29T 3 gTi - 3T 2 9 2 T 2 + 2gT 4 3T + 3 2 T 4 T 
-r- 1 (3T 9 2 T 3 - 2gT 3 9T - 4T + T 2 g 2 Tx - 6gT 2 gT! + 3g 2 T 2 T!) 

; _ r g2t dTo = 
or 

-at i + gg 2 T 2 - 2g 2 T 4 gti - 2gt 4 g 2 Ti + 2T 3 g 3 T 2 + 2t 2 g 3 T 3 (144) 
+r^ 1 (-2gti - 3gt 2 g 2 T 2 + t 2 g 3 t 2 + 2tig 3 t 3 - 2g 2 tigt 3 + g 2 f 4 gf + g 3 t 4 T ) 
: +r- 2 gf = o, 



6X A -A£ A = 

~A ~ . 

5ui — 5lo — 

6U-Au) = 



(143) 
1 



Or 
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(145) 

g 2 x 2 - gg 3 t 3 - 2gT 3 g 3 T 3 - 2r 3 g 4 T 3 - r 2 g 4 T 4 + 2gf 4 g 3 T 2 + 3g 2 t 4 g 2 T 2 
+r- 1 (+2g 3 T 2 gx 3 + 3g 2 x 3 g 2 T 2 - T 2 g 4 r 3 - g 4 x 4 Ti - 2g 3 T 4 gT 4 - g 2 T 4 g 2 T 4 ) 
-r- 2 g 2 T 4 = o, 



-g 3 f 3 + gg 4 T 4 - 2gt 4 g 4 T 3 - 4g 2 T 4 g 3 x 3 + 4g 4 t 4 gx 3 + 2g 5 r 4 f 3 (we) 
+r~ 1 (+2g 3 T 3 - eg 2 t 3 g 3 f 3 + g 3 t 2 g 2 f 4 + 3g 2 f 2 g 3 t 4 + 3gT 2 g 4 T 4 + x 2 g 5 t 4 - 2gt 3 g 4 T 3 ) 

: +r~ 2 g 3 T 2 = 0. 

These equations have the following relationship to each other: the applica- 
tion of the operator r 2 d r to the lower of each of the three consecutive pairs, 
[using Eqs. (|551) - (|101l) ]. leads to the equality with the operator g applied to the 
upper member of the pair, e.g., r 2 d r applied to Eq. (|146l) equals g applied to 
Eq.(Hig). 

Furthermore these four evolutionary Bianchi Identities are closely related 
to all the other evolution equations by similar identities. For example, from 
Eq.flSEl), i.e., 

5U- AS = (£ + 7-7)0; - v, (147) 

we have 

-gx 3 + gggx 4 + 2T 3 g 3 T 4 - 2gr 4 g 2 T 3 (i48) 
: +r- 1 (2gt 3 - 2gt 3 g 2 t 3 + t 2 g 3 t 4 + gf 2 g 2 x 4 ) + r~ 2 gT 2 = o. 

The application of g to (fHS]) yields Eq. ((M|) . 

~5v - AJi = Jl 2 + £(7 + 7) - 25/3 , (149) 

or 

(150) 

-g 2 x 3 + g 2 ggf 4 + 2t 3 g 4 T 4 + 2g 2 t 4 - 2g 2 f 4 g 2 t 3 - 2gr 4 g 3 T 3 + 2gr 3 g 3 T 4 
+r" 1 (2g 2 T 3 - 2gf 3 g 3 T 3 - 2g 2 t 3 g 2 t 3 + 2gf 2 g 3 t 4 + x 2 g 4 r 4 + g 2 x 2 g 2 x 4 ) 

: = -r~ 2 g 2 T 2 . 

Finally, in turn, the application of g to (|150l) yields the Bianchi Identity 
(|146[) . The commutator on a spin- weight— s function W s , i.e., 

mW s - mW s = -2sW s (151) 
has been used several times. 
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Returning to the evolutionary Bianchi Identities - leaving the first one, (|143|) , 
unchanged - we can see that the last three are simplified by performing the 
angular integrations. They can each be rewritten as 



5\Bi = 0, (152) 
d 2 B 2 = 0, (153) 
d 3 B 3 = 0, (154) 



with 



(155) 

Bx = -T' 1 +ggT2-2T2-2gT4STi+2T 3 3 2 T2-2gT 3 gT2 + 2g 2 T 3 T2 

+r~ 1 (f 2 g 2 r 2 - 2f i - 2gt 2 gf 2 + 2tig 2 t 3 - 2gtigr 3 + t g 2 r 4 ) + r~ 2 T , 



(156) 

b 2 = -t' 2 + ggx 3 - 2T 3 + T 2 g 2 T 4 - 2gT 4 gT 2 - gf 3 gf 3 + 2T 3 g 2 x 3 
+r- 1 (tig 2 T 4 - 2gT 3 gT 2 + x 2 g 2 T 3 ) + r- 2 Ti, 



B 3 = -T 3 +ggT 4 +2T 3 g 2 T 4 -2gT 4 gT 3 +r- 1 (2T 3 -gT 3 gT 3 +T 2 g 2 T 4 )+r- 2 T2, 

(157) 

so that they immediately integrate to 



Bi 


- Ki, 


(158) 


B 2 


= K 2 , 


(159) 


B 3 


= K 3 , 


(160) 



with (K l7 K 2 , K 3 ), the kernels of the operators (g, g 2 , g 3 ), i.e., (dK 1 ,(5 2 K 2 , (5 3 K 3 ) 
= 0. 

The kernels, using Eq. (p?T|) . are given by 

K 3 = A(() + B(()P + C(()P 2 +j-\D(0 + E(()P) + r- 2 F((), (161) 
K 2 = -D(C)-E(C)P-2r- 1 F(c), (162) 
K x = 2F(C). (163) 

with six arbitrary functions of (C). All have angular (or wire) singularities. Since 
we wish to restrict ourselves to regular solutions, the Ki are set to zero here. 
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Our final evolution equations for the potentials then are B3 = £? 2 = B\ = 
B = or 



Tn - 59Ti - 2T 3 g 2 T! + 23T 3 gf 1 
29T 3 gTo - 4T - 





-r- 1 (3f S 2 f 3 
r9 2 f 



t 2 S 2 Ti - 



(164) 

+ 29f 4 3t + g 2 T 4 T 
63T 2 5Ti + 3g 2 T 2 f!) 



-To 

dr 



4TpT 

or 



0, 



-f i + ggr 2 - 2T 2 



2gf 4 gf 1 + 2T 3 g 2 f 2 



2gx 3 gT 2 + 2g 2 T 3 x 2 



L (-2Ti + T 2 g 2 x 2 - 2gx 2 gT 2 



2Tig 2 T 3 - 2gTigx 3 + T g 2 r 4 



(165) 



2 T = 0, 



-t' 2 + ggf 3 - 2T 3 + T 2 g 2 T 4 - 2gT 4 gx 2 - gx 3 gT 3 + 2T 3 g 2 T3(i66) 

+r- 1 (T 4 g 2 T 4 - 2gT 3 gT 2 + T 2 g 2 T 3 ) + r~ 2 Ti = 0, 



(167) 

-T 3 + ggT 4 + 2T 3 g 2 T 4 -2gT 4 gT 3 + r" 1 (2T 3 -gT 3 gT 3 + T 2 g 2 T 4 )+r" 2 T 2 = 0. 

Finally, using Eqs. (|100j) and (|101[) of the ladder relations in Eq. (|167p . we 
obtain our non-linear wave equation T 4 which carries all the information of the 
evolution of the individual T„; 

d r % - r- 2 ggT 4 - <9 r 2 T 4 - 4r- 1 a r f 4 (168) 

= -29t 4 a r gf 4 + 45 r T 4 g 2 T 4 - r 5 r gt 4 <9 r gT 4 + r <9 2 t 4 g 2 f 4 . 

Since, with the ladder relations, Eq. (|168p is equivalent to the previous four 
equations, one could use either for integration. 

5 Examples 

A small set of solutions can be found by make the starting ansatz 

T 4 = c 4 - cgr- 1 + i C2 r- 2 - i Cl r~ 3 + ^c r~ 4 . (169) 
This leads immediately to 
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T - c , (170) 
Ti = ci-cor" 1 , 

T 2 = c 2 - cir + -c r , 

T 3 = c 3 - c 2 r _1 + \cxr~" 1 - ^c r~ 3 , 
2 o 

the c 's being independent of r. Inserting T 4 into (|167[) , and equating powers 
of r , yields the evolution equations for the individual c 's: 

c - 35ci - 2c 3 g 2 ci + 28c 3 3ci - 3c 2 3 2 c 2 + 23c 4 3c + c 3 2 c 4 = 0, (171) 

-c[ + 33c 2 - 2c 2 - 23c 4 3c! + 2c 3 3 2 c 2 - 23c 3 3c 2 + 2c 2 3 2 c 3 = 0, 

— c 2 + 33c 3 - 2c 3 + c 2 3 2 c 4 — 23c43c 2 — 3c 3 3c 3 + 2c 3 3 2 c 3 = 0, 

-C3 + 33c 4 + 2c 3 3 2 c 4 - 23c 4 3c 3 = 0. 

We obtain a system easily integrated if we make the further ansatz, 

c 4 = 0. (172) 

The ifTTTjl reduce to 

c' Q = 33ci + 2c 3 3 2 ci - 23c 3 3ci + 3c 2 3 2 c 2 , (173) 
ci = 33c 2 -2c 2 + 2c 3 3 2 c 2 - 23c 3 3c 2 + 2c 2 3 2 c 3 , 
c 2 = 33c 3 - 2c 3 - 8c 3 3c 3 + 2c 3 3 2 c 3 , 

4 = 0. 

leading to 

c 3 = c°(C,C), (174) 
c 2 = C ^(C,C)+«(33c!]-2^-3 C °3^ + 2 c 03 2 C °) 

ci = c°(C,C)+ J dw(33c 2 -2c 2 + 2c§3 2 c 2 -23c!]3c 2 + 2c 2 3 2 c§), 
c = c°(C, C) + / du(m Cl + 2c§3 2 ci - 23c°3ci + 3c 2 3 2 c 2 ). 

These solutions have in general cubic u dependence. 

The special case of c 3 = c 2 = c\ = 0, leaves Co = Cq(C, C)> a time independent 
solution 
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T = eg, (175) 
1 1 — c r , 



i D M — 2 



Another class of solutions of (|171|) starts with the ansazt, cq = c\ = C2 — 
C3 = 0. The last of (|171|) . becomes 

gg C 4 = 97(p 4 a c (F^ 2 c 4 ))-o. (176) 



The solutions, all type N, given by 

c 4 = P 2 G(C) + P 2 J p-^FiOdC, (177) 
are all singular, and have, unfortunately, wire singularities. 

6 Discussion 

In the interests of honesty, we describe the thought perturbations that led to 
the present investigation. Several months ago we were working on aspects of the 
complexified full Einstein equations - in particular - on asymptotic shear-free 
null geodesic congruences, and simply noticed that we could, with great ease, 
integrate the radial Bianchi Identities via the introduction of the five potentials 
(described in the text) with their ladder relations. Though we initially had no 
interest in investigating the self-dual Einstein metrics, the apparent simplicity, 
coming from the use of the potentials, charmed us into going for the "entire 
thing". Briefly it even seemed possible that we might go the entire way and 
be able to construct explicit vacuum self-dual metrics directly from the initial 
data. It turned out that this hope was very much dashed: first of all by the 
complexity of the remaining equations - they involved long and hard calculations 
- and more important, we ended with the unpleasant non-linear wave equation 
not easily solvable. 

So in the end the question remains: though we have found some pretty 
results, what did we really accomplish. These results certainly clarify the struc- 
ture of the self-dual equations - but do the results have applications or are they 
of interest in other investigations? They appear to potentially have use for our 
own investigations - but we are not certain. 

There are several further issues to be mentioned. 

The first is a mild mystery. The data needed for solutions to our non-linear 
wave equation are two complex functions, one the value of T4 on an initial null 
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cone uq, i.e., an arbitrary function of (r, £, £) and the second, C4, (the asymptotic 
value of T4) a news-like function of (u, £, £) that drives the evolution. On the 
other hand, from the original formulation of the self-dual Einstein metrics, via 
the so-called 'good-cut' equation, one needs only one complex function as the 
needed data, namely the asymptotic Bondi shear a°(uB, C; C); with ub the Bondi 
time. It is probable that by the proper counting of the ranges of the arguments 
it can be shown that the two sets contain the same amount of data. 

A second is that there should be a geometric meaning to the potentials, T„. 
We have not attempted to investigate that. 

For whatever it might be worth, and it is suggestive that it is worth some- 
thing, we point out that studies 18 of the self-dual Yang-Mills equations in the 
spin-coefficient form show that, for any gauge group, the entire set of Yang- Mills 
equations can be encapsulated into a single nonlinear wave equation very similar 
to our non-linear equation, (|168[) . The basic variable, #2, is one of three ma- 
trix valued potentials, foj all connected to each other by a ladder relationship 
analogous to the one for GR. 

The Yang-Mills wave equation 19 is 

dr& 2 ~ d^2 -2^8^ -r^fo = [SfoAfo], (178) 

with ladder relations among the potentials 

Si = r 2 d r % 2 , (179) 
do = r 2 d r $ 1 =r 2 d r {r 2 d r $ 2 ). (180) 

The tetrad components of the self-dual Yang-Mills field are given by 

Xo = r- 3 So, (181) 
Xi = r- 2 Qdi, (182) 
X 2 = r- 1 ^, (183) 

in an almost perfect analogy with the self-dual gravitational case. 
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